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ABSTRACT 


Technique: We present a method to determine the pressure at which significant cloud opacity 
is present between 2 and 6 bars on Jupiter. We use: a) the strength of a Fraunhofer absorption 
line in a zone to determine the ratio of reflected sunlight to thermal emission, and b) pressure- 
broadened line profiles of deuterated methane (CH3D) at 4.66 ^m to determine the location of 
clouds. We use radiative transfer models to constrain the altitude region of both the solar and 
thermal components of Jupiter’s 5-/rm spectrum. Results: For nearly all latitudes on Jupiter the 
thermal component is large enough to constrain the deep cloud structure even when upper clouds 
are present. We find that Hot Spots, belts, and high latitudes have broader line profiles than do 
zones. Radiative transfer models show that Hot Spots in the North and South Equatorial Belts 
(NEB, SEE) typically do not have opaque clouds at pressures greater than 2 bars. The South 
Tropical Zone (STZ) at 32°S has an opaque cloud top between 4 and 5 bars. From thermochemical 
models this must be a water cloud. We measured the variation of the equivalent width of 
CH3D with latitude for comparison with Jupiter’s belt-zone structure. We also constrained 
the vertical profile of H2O in an SEB Hot Spot and in the STZ. The Hot Spot is very dry for 
P<4.5 bars and then follows the H2O profile observed by the Galileo Probe. The STZ has a 
saturated H2O profile above its cloud top between 4 and 5 bars. 

Subject headings: planets and satellites: individual (Jupiter) — planets and satellites: atmospheres 


1. INTRODUCTION 


layers. Water was first detected on Jupiter by 


The abundance of H2O in Jupiter’s atmosphere 
is of fundamental importance in understanding the 
origin of Jupiter, the composition of its clouds, 
and Jovian dynamics beneath the upper cloud 



alyzed both KAO and Voyager IRIS spectra in 
Jupiter’s 5-/im spectral window to the deep at- 
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mosphere, finding a deep water abundance highly 
depleted with respect to the solar abundance. 
The spectrum between 4.5 and 5.4 provides 
a wealth of information about the gas composi¬ 
tion and cloud structure of the troposphere of this 
giant planet. Jupiter’s 5-/j,m spectrum is a mix¬ 
ture of scattered sunlight and thermal emission 
that changes significantly between belts and zones. 
Jupiter exhibits remarkable spatial structure at 
5 /im. Chemical models of Jupiter’s cloud struc¬ 
ture predict three distinct layers: an NH 3 ice cloud 
near 0.5 bars, an NH 4 SH cloud formed from a re¬ 
action of NH 3 and H 2 S at 2 bars, and a massive 
water ice/liquid solution cloud near 5 or 6 bars, de¬ 
pending on assu mptions of composition and ther - 
mal structure (seelWeidenschilling &: Lewid ( 1973h 
and IWong et al. ( 2015h '). Thermal emission from 
the deep atmosphere is attenuated by the variable 
opacity of one or more of these three cloud lay¬ 
ers. Hot Spots, located primarily in the North and 
South Equatorial Belts (NEB, SEB), exhibit 5- 
/rm radiances up to 70 times larger than surround¬ 
ing regions due to a minimum of cloud opacity. 
They also appear brighter than their surroundings 
at microwa ve wavelengths du e to a low ammonia 
abundance ( Sault et 111112004 ). 

The interpretation of 5-/im spectra of Jupiter, 
and the H 2 O abundance in particular, has been 
hampered by uncertainties in the pressure level of 
the lower boundary of the emitting region. Two 


differen t models have been proposed. iBioraker et al 


(Il986alli3l suggested that thermal radiation at 
5 ^m originates from levels as deep as 8 bars, 
310 K w here unit optic a l dep th in H 2 occurs. In 
contrast, ICarlson et al.l (1199211 proposed a model 
in which a massive water-ice cloud establishes 
the lower boundary near 5 bars, 273 K. Bjoraker 
et al. fitted KAO and Voyager IRIS spectra of 
H 2 O in Jupiter’s Hot Spots with a small abun¬ 
dance (4-30 ppm) distributed along a long path 
(60 km) between 2 and 8 bars. They also mea¬ 
sured H 2 O abundances in the low-flux zone regions 
using Voyager spectra. The NEB hot spots were 
found to be depleted in H 2 O between 2 and 4 bars, 
but belts, zones, and hot spots could all be fitted 
by the same H 2 O profile (4-30 ppm) betwee n 4 
and 8 bars. In marked contrast, Carlson et al.l fit¬ 
ted Voyager Hot-Spot spectra using a much larger 
mixing ratio (up to 3000 ppm, equivalent to 3 x 
the solar 0/H measured bv lAsolund et al.l ( 2009ll l 


confined to a narrow layer (10 km) between 4 and 
5 bars immediately above an opaque water cloud. 
Carlson et al. also examined Voyager spectra of 
the Equatorial Zone (EQZ) and other regions away 
from hot spots. All regions required a saturated 
H 2 O profile that increased from 300 ppm at 4 
bars to 3000 ppm at 5 bars, although NEB hot 
spots were sub-saturated in the 2 to 4 bar region. 


The Galileo Probe measured water vapor in 
situ at a single location on Jupiter, but unan¬ 
swered questions remain about its global abun¬ 
dance. The probe entered Jupiter’s atmosphere 
at 6.5°N planetocentric latitude near the south¬ 
ern portion of an NEB h ot spot. The Galileo 
Probe Mass Spectrometer (|Niemann et al.l Il992r) 
found water increasing with depth, from an up¬ 
per limit of 0.8 ppm at 2.7 bar, to measured 
mole fractions of 40±13 pp m at 11.0-11.7 bar and 
420: Jil40 ppm at 18-21 bar ( Niemann et al. (199^ 
and Wong et al. ( 2004bll l. The deepest value cor¬ 
responds to 0.45x solar 0/H. Ground-based imag¬ 
ing has shown that 5-^m hot spots, such as the 
one entered by the Galileo probe, co ver less than 
1% o f the surface area of Jupiter ( Orton et al.l 
1996ll . Thus, the H 2 O abundances observed by 


the probe may be characteristic of all or most hot 
spots, but they are probably not representative 
of Jupiter as a whole, especially since a range of 
other indirect studies (lightning flash depths, the 
tropospheric GO abundance, and discrete clouds 
at P > 4 ba r) point to solar o r supersolar water 
abundances (IWong et al.l 1200811 . Since the Galileo 
probe found that carbon , sulfur, and nitroge n were 
enriched by ^4x solar ( Wong et al. 2004bll . oxy¬ 
gen may also be enhanced by the same amount. 
Some formation models require ^lOx solar 0/H 
in order to trap Jupit er’s volatiles inside c ages of 
water ice clathrates ( Hersant et al.l l2004[) . The 
Juno mission, scheduled to begin orbiting Jupiter 
in July 2016, should answer many of these ques¬ 
tions. The Microwave Radiometer will measure 
water vapor bel ow Jupiter’s clouds to determine 
the 0/H ratio ( Janssen et al. 200511 . Interpreta¬ 
tion of these data may not be straightforward, 
however, due to e.g., the small microwave ab- 
sorptivity of HgO g as compared with NH 3 (see 


de Pater et al.l (1200511 for details). A complemen¬ 


tary way to determine the deep H 2 O abundance, 
therefore, is highly desirable. The Jovian In¬ 
frared Auroral Mapper (JIRAM) will acquire near- 
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infrared spectra of Jupiter, in cluding the 5-micron 


window ( Adriani et al. 20l4l . 


Water ice has been detected in isolated regions 
where active convectio n lofted the ice well abov e 
its condensation level ( Simon-Miller et al. 2000ll . 


but water clouds are generally h idden by overlying 
NHfi and NH 4 SH clouds (e.g. ISromovskv fc Fry 
(l2ninll l. Only very thin cl ouds were found i n 


1998I1 . 


the Galileo Probe Hot Spot (|Ragent et al 
consistent with condensation via weak turbu¬ 
lent updrafts within the descending b ranch of 
an eq u atorially-trap p ed R ossby wave (jFriedson 


Wong et al. ( 201,^). M odels of NIMS 
by Nixon et al.l ( 200ll l in clude water 


(120051), 
spectra 

clouds in at least some Hot Spots, and iRoos-Serote et al.l 
(2004) showed that NIMS spectra of Hot Spots 
cannot rule out water clouds whose opacity is 
entirely restricted to P > 5 bar. Evidence for 
wate r clouds in Jupiter’s zones is also ambigu¬ 
ous. IProssart et al. ( 1998 ) compared dayside and 
nightside NIMS spectra of low-flux regions in the 
EQZ. A saturated H 2 O profile above an opaque 
water cloud at 5 bars provides a satisfactory fit 
to these data. At the spectral resolution of IRIS 
(4.3 cm~^) and NIMS (10 cm“^) we can only 
retrieve a column abundance of H 2 O above an 
assumed lower boundary, which can be either a 
water cloud or opacity due to H 2 . Thus, we sim¬ 
ply cannot tell whether or not Hot Spots or zones 
have water clouds. This problem will also affect 
the interpretation of 5-/im spectra from JIRAM 
on Juno, which has the same s pectral resolution 
as Voyager/IRIS (|Adriani et al. 


B Pectr al 

lioli. 


In Section 2 we present ground-based observa¬ 
tions of Jupiter’s belts and zones that have suffi¬ 
cient spectral resolution to resolve line shapes. We 
demonstrate that we can derive the cloud struc¬ 
ture at P > 2 bars even when higher-altitude 
clouds greatly attenuate the thermal flux from 
the deep atmosphere. We used line shapes to 
detect water clouds and to determine the pres¬ 
sure at which these clouds become optically thick. 
This resolves the ambiguity of whether the lower 
boundary of the 5-/im line formation region is due 
to H 2 or due to opaque clouds. This, in turn, will 
yield more accurate gas abundances for use in con¬ 
straining models of Jupiter’s origin and in under¬ 
standing the dynamics of the atmosphere beneath 
the upper clouds. 


2. OBSERVATIONS 


Five-micron spectra of Jupiter were acquired 
using NIRSPEC on the Keck 2 telescope on March 
11, 2014. NIRSPEC is an echelle spectrograph 
with 3 orders dispersed onto a 1024x1024 InSb ar¬ 
ray at 5 /im at our selecte d grating/cross-disp erser 


settings of 60.48 / 36.9 ( McLean et al. 19981) . A 


0.4" X 24"slit was aligned north-south on the cen¬ 
tral meridian of Jupiter, resulting in spectra with 
a resolving power of 20,000. NIRSPEC has an ad¬ 
vantage over instrumentation on other telescopes 
(for example CSHELL on the IRTF) due to the 
fact that 3 echelle orders at 5 ^m are placed on 
the detector array. Thus, each pixel along the slit 
corresponding to different latitudes on Jupiter has 
simultaneous spectra at 4.6, 5.0, and 5.3 /rm. In 
this paper we will focus on Order 16 covering 2131- 
2165 cm“^ (4.62-4.69 /rm) and Order 15 which 
covers 1999-2031 cm“^ (4.92-5.00 /im). Jupiter 
subtended 41" and the geocentric Doppler shift 
was 26.3 km/sec. The water vapor column above 
Mauna Kea was 2 precipitable mm derived from 
fitting telluric lines in the Jupiter spectra. We did 
not use stellar spectra for flux calibration or at¬ 
mospheric transmission because the humidity dou¬ 
bled between the time of the Jupiter and stellar 
observations. The flux calibration is described in 
Section 3. Fig. 1 shows an image of Jupiter us¬ 
ing the SCAM guide camera on NIRSPEC. Al¬ 
though the spectroscopy was performed at 5-/im, 
the guide camera works at shorter wavelengths. A 
K-prime filter centered at 2.12 /rm was used to 
obtain sufficient contrast to separate the bands of 
variable haze reflectivity overlying Jupiter’s belts 
and zones. Two slit positions were required to ob¬ 
tain pole-pole spectra of Jupiter. In Eig. 1 spatial 
pixels that exhibit maxima in flux at 4.66 /rm are 
shown in red, locations that exhibit narrow line 
profiles at 4.66 /im are shown in blue, and regions 
that have broad pressure-broadened line profiles 
are shown in green. We focus on two spatial loca¬ 
tions with characteristically different spectra: Re¬ 
gion A in the South Tropical Zone (STZ) at 32°S, 
and Region B, a Hot Spot in the SEB at 17°S. We 
developed radiative transfer models to calculate 
synthetic spectra for regions A and B, as described 
below. 

The band pass at 4.66 /im contains absorption 
lines of deuterated methane (CH 3 D), phosphine 
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Fig. 1.— Keck/SCAM image of Jupiter at 2.12 fim shows 
2 NIRSPEC slit positions, one in the northern, and one in 
the southern hemisphere. Curvature is due to navigating 
positions onto an image taken a few minutes earlier. Red 
pixels denote Hot Spots. Blue pixels denote positions of 
minima in deuterated methane (CH 3 D) equivalent width 
and thus candidates for opaque clouds. Green regions de¬ 
note maxima in CH 3 D equivalent width and thus regions 
without opaque clouds. Radiative transfer models were 
used to calculate synthetic spectra for the zone marked A 
and the SEB Hot Spot labeled B. 


(PH 3 ), and H 2 O. Methane and its isotopologues 
do not condense, and they are not destroyed pho- 
tochemically, in Jupiter’s troposphere. We there¬ 
fore assume that CH 4 and CH 3 D have a con¬ 
stant mixing ratio with respect to H 2 in Jupiter’s 
troposphere, which means that variations in the 
strength and shape of CH 4 and CH 3 D lines be¬ 
tween belts and zones on Jupiter should be due to 
changes in cloud stru cture, not gas concentration. 


Bioraker et al.l (|2002ll measured the spatial varia¬ 


tion of the weak 1 / 3 - 1^4 band of CH 4 at 5.18 /rm in 
an attempt to derive cloud structure. Unfortu¬ 
nately, they were unable to calculate synthetic 
spectra that ht this feature. This is possibly due to 
inaccurate or incomplete spectroscopic parameters 
such as line strengths and broadening parameters. 
Spectroscopic parameters for the much stronger 
1^2 fundamental band of CH 3 D at 2200 cm~ ^, or 
4.5 /rm, are well known (INikitin et al.Nl997l) and 
as shown below we now are able to derive cloud 
structure and spatial variations therein. 

Fig. 2 shows spectra at 4.66 (2144 cm“^, 


Variation in CH^D Line Profiles 



Wavenumber 

Fig. 2.—^ Jupiter spectra at 4.66 pm (NIRSPEC order 
16) show dramatic variation with latitude in the strength 
and width of CH 3 D lines. Four absorption features are de¬ 
noted by vertical lines. Narrow CH 3 D lines are observed in 
the Equatorial and North Tropical Zones (EQZ, NTZ), but 
they are much broader in Hot Spots and at high northern 
latitudes. A phosphine (PH 3 ) line and a Doppler-shifted 
H 2 O line are shown. Fraunhofer lines (CO in the Sun) are 
only observed in low-flux zones. T denotes telluric lines. 


NIRSPEC order 16) of 5 regions in Jupiter’s north¬ 
ern hemisphere that exhibit either a maximum 
in flux (an NEB Hot Spot at 8.5°N), minima 
in the equivalent width of CH 3 D (the EQZ at 
0.5°N and NTZ at 23°N), or local maxima in 
CH 3 D equivalent width. A horizontal bar indi¬ 
cates the limits that we chose for numerical in¬ 
tegration of the equivalent width that includes 
six absorption lines of CH 3 D (blended into four 
at this spectral resolution). All spectra are nor¬ 
malized to 1.0 at 2141.6 cm“^ to facilitate com¬ 
parison of line shapes, since the radiance of the 
NEB Hot Spot is 70 times that of the NTZ. A tel¬ 
luric H 2 O line, its Doppler-shifted counterpart, 
a PH 3 feature, and a Fraunhofer line due to CO 
in the Sun are also shown. Note that the Fraun¬ 
hofer line is observed only in low-flux zone regions 
such as the EQZ and NTZ. In Sections 3 and 4 
we describe how the strength of this feature can 
be used to constrain cloud models. The equivalent 
width of the set of CH 3 D lines at 2144 cm“^ varies 
dramatically from belt to zone. Hot Spots and 
belts exhibit broad pressure-broadened line pro¬ 
files, while zones have much narrower features 
that are resolved as 4 distinct CH 3 D absorption 
lines. Note that an isolated CH 3 D line would be 
spectrally resolved at a resolving power of 20,000 
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(0.11 cm“^ resolution). This is because molecules 
such as CH 3 D typically have broadening coeffi¬ 
cients of 0.06 cm“^/atm and the line formation 
region on Jupiter at 4.66 /rm takes place at pres¬ 
sures greater than 2 bars. 


3. DATA ANALYSIS AND RESULTS 


Absolute flux calibration was performed for 
both NIRSPEC orders 15 and 16 by smoothing 
the spectrum of the NEB Hot Spot at 8.5°N to 
4.3 cm~^ resolution, dividing by the transmittance 
of the Earth’s atmosphere above Mauna Kea, and 
scaling the radiance of the resulting spectrum to 
an average of the 4 hottest spectra of Jupiter’s 
NEB observed by Voyager IRIS in 1979. 


In addition to pressure-broadened line pro¬ 
files, the spectrum at 4.66 /rm provides valu¬ 
able information on the ratio of reflected sun¬ 
light to thermal emission on Jupiter, which is 
critical for understanding Jupiter’s cloud struc¬ 
ture. We compared the equivalent width of the 
Eraunhofer line at 2141.8 cm“^ with its measured 
value in the Sun using data from the Atmospheric 
Trace Molecule Spectroscopy (ATMOS) experi¬ 
ment whi ch flew on the space shuttle Challenge r 
in 1985 ( Earmer fc Norton Il989t iFarmed Il994 ). 
Solar Eraunhofer lines in spectra of Jupiter arise 
from regions with a significant fraction of reflected 
sunlight, such as locations with thick high clouds. 
Eraunhofer lines are not observed in regions that 
are dominated by thermal emission, such as Hot 
Spots. We selected a zone at 32°S (marked“A” 
on Figures 1, 4, and 5) for further study based 
on the strength of this particular Fraunhofer line. 
It is 43% as strong as in the Sun; thus, 57% of 
the flux consists of thermal emission originating 
in the deep atmosphere that has been attenu¬ 
ated by one or more cloud layers before escaping 
to space. This thermal flux preserves the broad 
line profiles of CH 3 D caused by collisions with 
H 2 and He in its line formation region at the 
deepest levels probed. We excluded the spectrum 
of the NTZ shown in Fig. 2 for studies of the deep 
atmosphere because the Fraunhofer line is 95% 
as strong as in the Sun; thus, there is insufficient 
thermal flux to constrain the deep cloud structure 
at 23°N. We also modeled the SEB Hot Spot at 
17°S and 208°W (marked ”B” on Figures 1, 4, and 
5), which had the highest radiance in our entire 


dataset. The noise level was measured using the 
standard deviation of the number of counts in the 
spectral pixel corresponding to a saturated telluric 
water line in each order. This was evaluated over 
11 spatial pixels over the lowest flux region (the 
STZ) and also for 11 spatial pixels off of Jupiter. 
Both gave the same result. In Order 16 the sig¬ 
nal to noise ratio (S/N) of the SEB Hot Spot at 
2141.6 cm“^ was 2500; the S/N of the zone at 
32°S was 48. In Order 15 where strong H 2 O lines 
occur, the S/N of the same Hot Spot was 1900 at 
2012.2 cm-i ; the S/N of the STZ was 35. 

Synthetic spectra were calculated using the 
Spectrum Synthesis Progra m (SSP) radiative 
transfer code as described in Kunde fc Maguir j 


1119741) . The input temperature profile was ob¬ 
tained from the Galileo Probe (Seiff et al. 1998). 
Line parameters for CH 3 D and other 5-um ab¬ 
sorbe rs are from GEISA 2003 (jjacauinet-Husson et al 
2005 ). Parameters for CH 3 D-H 2 and CH 3 D- 
He broadening have been measured in the lab 


He broadening nave been mea,sureQ m tne lab 
( Boussin et al.lllO^ Lerot et al.l2003t Feiard et al. 


20031 ). We used a broadening coefficient of 
0.0613 cm-i /atm (296/T)0'5 for CH 3 D colliding 
with a mixture of 86.3% H 2 and 13.6% helium, as 
measured by the Galileo Probe (von Zahn et al. 
1998). Pressure-induced H 2 coefficients were ob- 
tained using laborat ory measurements at 5 /rm by 
lBachet_e£jh J 1983[) an d the f ormalism developed 
by Birnbaum fc Goh^ ( 1976[) . 


The base of the model was set to 20 bars, 416K 
for the SEB Hot Spot to ensure that the base was 
well below the level where Jupiter’s atmosphere 
becomes optically thick due to pressure-induced 
H 2 opacity. For the zone at 32°S, we investigated 
lower boundaries at 2, 4, 5, and 20 bars to simu¬ 
late opaque NH 4 SH clouds, H 2 O clouds, and no 
deep clouds in the troposphere. For the Hot Spot, 
we calculated a spectrum free of deep clouds for an 
emission angle of 18.7° andforanH 2 0 mole frac¬ 
tion of 47 ppm for P > 4.5 bars as shown in Fig. 3, 
consistent with res ults from the Galileo Probe 
( Wong et al.l 12004^ . The model was iterated to 
ht CH 3 D (0.18 ppm) and PH 3 (0.45 ppm). Our 
ht to CH 3 D in the Hot Spot is ve ry close to the 
value (0.16±0.04 ppm) derived bv iLellouch et al. 
(I 2 OOII) using spectra at 8.6 /rm from the In¬ 
frared Space Observatory. We then convolved 
the Jupiter spectrum to 0.02 cm“^ , doppler- 
shifted it by 25.2 km/sec, multiplied by the calcu- 
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lated transmittance above Mauna Kea, convolved 
it to 0.14 cm“^ to match the NIRSPEC data, 
and multiplied the spectrum by 0.332 to simu¬ 
late the transmission of upper cloud layers on 
Jupiter to match the observed continuum. As 
shown in Fig. 4, the synthetic spectrum fits Re¬ 
gion B, the SEB Hot Spot spectrum without any 
deep cloud opacity (P > 2 bars). All 4 CH 3 D fea¬ 
tures, including the wing between 2143.2 and 
2143.6 cm“^ are matched. 

The zone model is more complicated. A ra¬ 
diative transfer model was used to calculate spec¬ 
tra for an emission angle of 33.4°. Above each 
lower boundary, the gas composition was set to 
0.18 ppm CH 3 D. For the zone model with no 
deep clouds, the Galileo Probe H 2 O profile was 
used. For models with opaque clouds at 2, 4, and 
5 bars, a saturated profile of H 2 O was used (see 
Fig. 3). The mole fraction of PH 3 was iterated to 
a value of 0.7 ppm to match the absorption fea¬ 
ture at 2143 cm“^ for lower boundaries at 4 and 
5 bars. For the zone model without deep clouds, 
the PH 3 abundance in the STZ was iterated to 
a value of 0.45 ppm. Thus, PH 3 can also be 
used to discriminate between deep cloud models, 
as we describe below. A larger PH 3 mole frac¬ 
tion (0.7 ppm) above an opaque cloud layer at 4 
bars yields the same PH 3 absorption as a smaller 
mole fraction (0.45 ppm) above the (deeper) level 
where the atmosphere becomes opaque due to H 2 - 
H 2 opacity. The synthetic spectrum was split into 
two parts. The reflected solar component was cal¬ 
culated using the transmittances above a reflect¬ 
ing layer at 300 mbar and for an upper cloud re¬ 
flectance of 9%. It was convolved to 0.02 cm“^ res¬ 
olution, doppler-shifted, and multiplied by the 
ATMOS solar spectrum. The thermal component 
was convolved, doppler-shifted, and multiplied by 
transmittances of 0.393, 0.0269, 0.016, and 0.0038 
for models with a base of 2, 4, 5, and 20 bars, re¬ 
spectively, and added to the reflected spectrum. 
This, in turn, was multiplied by the transmis¬ 
sion above Mauna Kea and finally convolved to 
0.14 cm“^ resolution. 

Due to the numerous parameters in this model, 
we illustrate the parameter sensitivity of the 
model for the reflected solar and thermal com¬ 
ponents of the zone model separately in Fig. 5. 
We treat the reflective upper cloud as spectrally 
grey. Ice components such as NH 3 , NH 4 SH, 


and H 9 O have bee n seen at other wavelengths 


(e.g., 

Brooke et al. 

(19981. Simon-Miller et al. 

(29991. 

Baines et al.l ( 

299211. Wong- et al.l (2994a 1. 


Sromovskv fc FrvI (I 2 OIO II. However, our 


spec¬ 
tral windows do not include significant ice ab¬ 
sorption features. The imaginary indices of re¬ 
fraction o f NH3 ( Marto nchik et Il984l l and 


299711 are on the order 


NH 4 SH (|Howett et al. 

of 19“^ at the wavelengths studied here, much 
lower than values closer to unity at the 1^3 vibra¬ 
tion transitions for these ices (at 3.9 and 3.4 /im, 
respectively). Thus, the spectral features in the 
reflected solar component are due to CH 3 D ab¬ 
sorption lines in the upper troposphere of Jupiter 
or Fraunhofer lines in the Sun. 


In the top panel of Fig. 5 we show 4 Fraunhofer 
lines, marked S. The strength of the strongest Jo¬ 
vian feature at 2141.8 cm“^ when compared with 
the spectrum of the Sun acquired by the ATMOS 
investigation yields the relative fractions of re¬ 
flected sunlight (9.43) and thermal emission (9.57) 
for this portion of Jupiter’s South Tropical Zone. 
The fraction of reflected sunlight combined with 
the calibrated continuum level allows us to derive 
an upper cloud reflectance of 9%. There are 5 
CH 3 D absorption features in this spectral range. 
Ammonia clouds and hazes limit the penetration 
of reflected sunlight to the upper troposphere. We 
investigated the pressure of the reflecting layer by 
calculating spectra at 199, 399, and 699 mbars. 
We adopted a value of 390 mbars for the STZ. 
Note that in order to match the narrow absorption 
cores of CH 3 D in a model that consists of the sum 
of two components, there are a family of solutions 
for the pressure level of the solar reflecting layer 
and for the thermal cloud top pressure. Increasing 
the pressure of the reflecting layer from 300 to 600 
mbars may be compensated by, for example, de¬ 
creasing the (deep) cloud top pressure from 4 bars 
to 2 bars in the thermal component. However, we 
can exclude this possibility by studying the wing 
of the CH 3 D feature between 2143.2 and 2143.6 
cm“^. The top panel shows that the reflected solar 
spectrum is flat over this range. In contrast, each 
thermal model has a different CH 3 D wing line 
slope. We now return to Fig. 4 which compares 
the sum of the reflected and thermal components 
of each model to the observed STZ spectrum. We 
can exclude models with an opaque NH 4 SH cloud 
at 2 bars (blue curve) as well as the model with no 
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deep cloud (red curve). Note that the calculated 
radiances for both of these models lie well outside 
the error bars of the zone spectrum. The model 
with a deep cloud at 5 bars fits portions of the 
spectrum but the model with an opaque cloud at 
4 bars (green curve) provides a better overall fit 
to the spectrum. 

The NIRSPEC spectral bandpass for Order 16 
covers more than just this particular CH 3 D ab¬ 
sorption feature. In Fig. 6 we compare the STZ 
spectrum with the same zone models for a spec¬ 
tral region adjacent to the one displayed in Figs. 
4 and 5. Here there are three telluric H 2 O lines. 
To the left of each telluric line is a Jovian water 
line red-shifted by 0.18 cm“^ due to the relative 
velocity of Jupiter with respect to the Earth. The 
wings of these water lines permit us to discrim¬ 
inate between models. The best fit requires an 
opaque cloud between 4 and 5 bars. 


Next, we compare our retrieved PH 3 abun¬ 
dances in the 4 to 8 -bar level of the STZ with 
measurements of PH 3 at 1 bar from the Cassini 
flyby of Jupiter in January 2001. lirwin et al. 
(|2004ll retrieved PH 3 mole fractions on Jupiter 
ranging from 0.9 to 1.5 ppm using zonal aver¬ 
ages of CIRS spectra at 9 /rm between 60°S and 
60°N. The retrieved PH 3 mole fraction at 32°S 
was 1.0±0.2 ppm. Phosphine falls off with height 
in the upper troposphere due to ultraviolet pho¬ 
tolysis, but its mole fraction is not expected to 
change between 1 and 8 bars. There remain dis- 
crepa ncies between PH 3 retrievals at 5 and 9 


/rm (jFletcher et al.l I 2 OO 9 II . Nevertheless, better 
agreement for the abundance of PH 3 at 32°S is 
achieved for a model with opaque clouds near 4 to 
5 bars (0.7 ppm) than for a zone model with no 
deep cloud (0.45 ppm). 


Thus, using three independent arguments that 
are based on: a) the slope of the CH 3 D line wings, 
b) the slope of the H 2 O line wings, and c) the 
derived mole fractions of PH 3 , we conclude that 
there must be significant cloud opacity between 4 
and 5 bars in the STZ at 32° S. Based on the tem¬ 
peratures at these pressure levels (257K to 275K), 
thermochemical models predict that this is a water 
cloud. 


In Fig. 7 we compare the spectrum of the SEB 
Hot Spot marked B at 2016 cm“^ (4.96 fim) with 3 
different models calculated using the vertical pro¬ 
files of H 2 O shown in Fig. 3. The parameters 


are the same as those used to fit the spectrum at 
2144 cm“^ except for H 2 O and NH 3 . We initially 
used a vertical prohle of NH 3 derived from ab- 
sorption of the radio signal from the Galileo Probe 
( Folkner et al. (199^ as modified bv iHanlev et al 


(120091) 1. We found that a scaling factor of 2 was 
required to ht the NH 3 absorption features shown 
in Fig. 7. This corresponds to NH 3 mole frac¬ 
tions ranging from 240 ppm at 2 bars to 600 ppm 
at 4.5 bars. This is significantly larger than val¬ 
ues derive d from ground-base d microwave observa¬ 
tions (see lSault et al.l ( 2004h l. However, we have 
not used the much stronger absorption features 
in Order 14 (5.3 /rm) to constrain the NH 3 ver¬ 
tical profile. Thus, the NH 3 mole fractions re¬ 
ported here should be regarded as preliminary. 
The three different vertical prohles of H 2 O in 
Fig. 3, and used in our calculations, were based on 
the mass spectromete r data on the Galileo Probe. 
Niemann et al. ( 1998ll report ed an upper limit to 


H 2 O of 0.8 ppm at 2.7 bars. Wong et al. ( 2004bh 


reported a mole fraction of H 2 O of 40±13 ppm at 
11.0-11.7 bars on Jupiter. The only adjustable pa¬ 
rameter in our model was the pressure at which the 
H 2 O mole fraction increased from 0.8 to 47 ppm, 
as shown in Fig. 3. The best ht is for a pressure 
of 4.5 bars. Assuming that the SEB Hot Spot is 
similar to the one that Galileo entered, this data 
point provides a useful measurement of the depth 
at which dynamical processes have dried out Hot 
Spots on Jupiter. This depth compares w ell with 
the depth determined bv ISault et al.l ( 2004h for the 
NH 3 abundance in hot spots from microwave ob¬ 
servations. 


We next modeled the zone marked A at 
2016 cm“^. We used the same procedure as was 
used to model the zone spectrum at 2144 cm“^. 
We explored the same set of models with cloud 
tops at 2, 4, and 5 bars, as well as the model with 
no deep cloud. This portion of Jupiter’s spectrum 
is not sensitive to the deep cloud structure. Using 
the results from htting CH 3 D, we adopted a model 
with an opaque cloud top at 4.0 bars and a satu¬ 
rated H 2 O prohle above. The NH 3 mole fraction 
was iterated until a rough ht was achieved using 
a value of 250 ppm between 0.75 and 4.0 bars and 
a saturated value for P < 0.75 bars. Without 
any obvious Fraunhofer lines in this spectral re¬ 
gion to constrain the rehected solar component, 
we assumed the same upper cloud rehectance of 
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9% as at 2144 cin“^ and we obtained an upper 
cloud transmittance of 0.0542. The transmittance 
is therefore less certain than the value derived at 
2144 cm“^. The principal conclusion from this 
spectral region is that the observed H 2 O line 
profiles in the STZ are slightly broader than the 
telluric H 2 O features, but narrower than H 2 O fea¬ 
tures in the SEB Hot Spot. A saturated H 2 O pro¬ 
file for P < 4.0 bars provides a satisfactory fit to 
the observed spectrum. 



Fig. 3.— Vertical profiles of H 2 O used to calculate syn¬ 
thetic spectra in Figures 4 and 5. 










2.5 10 '^ 


vio 


Fig. 4. — The pressure level of deep clouds on Jupiter as constrained by CH 3 D line profiles observed in Order 16. An SEB 
Hot Spot (gray curve) at 17°S (labeled B in Fig. 1) is fitted without any clouds at pressures greater than 2 bars (pink curve). 
Hot Spot error bars would be smaller than the thickness of the gray curve. The spectrum of the STZ at 32° S (points with 
error bars, and labeled A in Fig. 1) was modeled using opaque clouds at 2, 4, 5 bars and a model without any deep clouds. We 
may exclude an opaque NH 4 SH cloud at 2 bars due to the poor fit of the blue curve to the observed spectrum. The best fit 
requires an opaque cloud between 4 and 5 bars. Thermochemical models predict that cloud opacity at this pressure level is due 
to a water cloud. Note the factor of 50 difference in radiance scales between the Hot Spot and the zone spectra. The calculated 
transmittance above Mauna Kea for 2 mm precip H 2 O is shown as a dashed blue line. T denotes a telluric H 2 O line and S 
denotes a solar (Fraunhofer) line. 
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Fig. 5. — Top: The reflected solar component of a radiative transfer model of the STZ at 32°S. Four Fraunhofer lines are 
marked as S. A telluric H 2 O line is marked as T. A reflecting layer is placed at 100, 300, and 600 mbars in Jupiter’s upper 
troposphere. One PH 3 and five CH 3 D lines are shown. The narrow cores of the CH 3 D lines are sensitive to the solar 
component. Bottom: The thermal component is shown for models with cloud tops at 2, 4, and 5 bars, and for a model with 
no deep cloud. The wing of CH 3 D between 2143.2 and 2143.6 cm~^ and the PH 3 line at 2143.0 cm~^ are due to the thermal 
component. Radiance scales are offset to match the continuum of the observed spectrum. 
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Fig. 6 .— The spectrum of the STZ at 32°S and radiative transfer models with opaque clouds at 2, 4, and 5 bars, and a model 
without any deep clouds are shown for a spectral region adjacent to the CH 3 D features. T denotes telluric H 2 O lines and 
S denotes a solar (Fraunhofer) line. Note the 3 Jovian H 2 O lines red-shifted by 0.18 cm~^ from their telluric counterparts. 
The wings of these water lines permit us to discriminate between models. The best fit requires an opaque cloud between 4 and 
5 bars. 
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Fig. 7. — Vertical profiles of water vapor mixing ratio as constrained by H 2 O line profiles observed in Order 15. An SEB Hot 
Spot (Region B, gray curve) at 17°S is fitted without any deep clouds. Hot Spot error bars would be smaller than the thickness 
of the gray curve. Three vertical profiles of H 2 O were calculated for the Hot Spot. The best fit used the Galileo Probe value 
for H 2 O for pressures greater than 4.5 bars (see also Fig. 3). The spectrum of the STZ at 32°S (Region A, points with error 
bars) was fitted using an opaque cloud at 4 bars and a saturated H 2 O profile above it (purple curve). The H 2 O lines in the 
zone are slightly broader than the telluric features, but much narrower than in the Hot Spot. 
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Fig. 8 .— The flux at 4.66 /im (red curve), the equivalent width of CH 3 D (black curve), and the the ratio of the equivalent 
width of the Fraunhofer line at 2141.8 cm~^ on Jupiter to that in the Sun is shown for each latitude along the central meridian 
on Jupiter. Hot Spots are present at 8.5°N and 17°S. Latitudes exhibiting minima in CH 3 D equivalent widths are candidates 
for water clouds. Radiative transfer models were used to match a zone spectrum at 32°S (marked ”A”) and a Hot Spot spectrum 
at 17°S (marked ”B”). The gold curve indicates the fraction due to reflected sunlight. Large values imply thick upper clouds 
that attenuate thermal radiation from below and reflect sunlight above. The bar on the right denotes haze reflectivity variation 
with latitude obtained from 2.12 /im SCAM images, roughly similar to belt-zone structure seen at visible wavelengths. 
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4. DISCUSSION 

This paper describes a technique to determine 
water cloud heights using spectrally-resolved line 
shapes of CH 3 D. We have demonstrated with a set 
of spectra from characteristic spatial regions that 
variation in the water cloud height manifests as a 
variation in the equivalent width of CH 3 D. This 
technique can be applied to all regions covered by 
the spectrometer slit. We show the CH 3 D equiv¬ 
alent width as a function of latitude (black curve) 
in Fig. 8, highlighting the latitude regions which 
may have water clouds. The observations are ex¬ 
tremely valuable both for understanding Jupiter’s 
atmospheric circulation, as well as its bulk water 
abundance (a constraint on planetary formation 
scenarios). 

However, the equivalent width of CH 3 D alone 
is not sufficient. In order to map H 2 O clouds 
on Jupiter, we need to measure the latitudinal 
profile of both CH 3 D line shapes as well as the 
strength of the Fraunhofer line at 2141.8 cm“^. 
In Fig. 8 we plot three quantities as a function 
of latitude on Jupiter. First is the radiance inte¬ 
grated over the 4.6-/im bandpass and normalized 
to 1.0 as a function of spatial pixel and converted 
to planetographic latitude. This is shown in red. 
The next quantity is the equivalent width of our 
set of CH 3 D absorption lines integrated between 
2143.21 and 2145.04 cm“^ (denoted by the hori¬ 
zontal bar in Fig. 2). This is shown in black. The 
last quantity plotted in Fig. 8 is the fraction of 
reflected sunlight as a function of latitude. This 
was obtained by measuring the equivalent width 
of the Fraunhofer line at 2141.8 cm“^ as a func¬ 
tion of latitude and dividing it by the value in the 
Sun as measured by ATMOS. At latitudes such as 
23°N where this value is close to 1.0, there is in¬ 
sufficient thermal flux to constrain the deep cloud 
structure. However, by using adjacent latitudes 
where CH 3 D equivalent widths are also small, and 
where the fraction of reflected solar flux < 0 . 6 , we 
can infer the presence of water clouds. Thus, we 
have indicated ’’possible water clouds” for the en¬ 
tire region where values of CH 3 D equivalent width 
are small. The vertical bar on the right denotes 
the brightness of Jupiter in reflected sunlight at 
2.12 ^m obtained from SCAM images navigated 
onto a planetographic latitude grid. This wave¬ 
length sounds hazes in the upper troposphere and 


illustrates Jupiter’s belt-zone structure. 

The equivalent width of CH 3 D is sensitive to 
two parameters: the fraction of reflected sunlight 
and the pressure level of the deep cloud. Consider 
two regions on Jupiter. They both lack an opaque 
water cloud, but one has thin upper clouds so that 
the 5-/xm spectrum is 100% thermal. The other 
has a thick, reflective upper cloud that attenuates 
the thermal component yielding a spectrum that 
is 50% thermal and 50% reflected solar. In the 
latter case, the contribution of the reflected solar 
continuum will significantly reduce the equivalent 
width of the CH 3 D absorption features that are 
present in the thermal component. A low value 
of CH 3 D equivalent width might be interpreted 
as due to a deep water cloud, whereas in this case 
the small equivalent width was due solely to the 
presence of the solar component from an upper 
reflective layer. Thus, we have to be cautious in 
using the equivalent width of CH 3 D alone to in¬ 
fer deep cloud structure. Nevertheless, the data 
plotted in Figure 8 show interesting structure at 
latitudes where the fraction of reflected solar radi¬ 
ation is essentially zero. High latitudes, southward 
of 55°S and northward of 35°N show interesting 
variations in CH 3 D equivalent widths that can¬ 
not be attributed to reflected sunlight and thus 
are related to deep cloud structure. 

With the caveat that radiative transfer mod¬ 
eling is required to verify the presence of wa¬ 
ter clouds in regions with small CH 3 D equiva¬ 
lent widths and significant fractions of reflected 
sunlight, we can make the following conclusions: 
Where thick water clouds are present at pres¬ 
sures in the 3-5 bar range, the CH 3 D equiva¬ 
lent width is small, and the black curve (Fig. 8 ) 
tends toward lower values. Regions with thick, 
high-altitude (P < 4 bars) water clouds can be 
associated with widespread upwelling, since down- 
welling would provide the opposite effect of clear¬ 
ing cloud opacity via sublimation. We find the 
smallest CH 3 D equivalent widths in the three low- 
latitude zones: the equatorial zone and the north 
and south tropical zones. This implies a large ver¬ 
tical extent to the upwelling responsible for form¬ 
ing the thick white clouds observed in zones at 
optical wavelengths. 

Thick water clouds in zones specifically con¬ 
tradict models of inverted two-la yer circulation 
within the tropospheric cloud decks, lingersoll et al. 
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(l 200 (Tl) and lShowman fc de PateJ (|2005l) suggested 
inverted two-layer circulation schemes to explain 
widespread ammonia gas depletion even in zones. 
In this type of circulation model, regions of up- 
welling (thick clouds) in the visible upper levels 
(NH 3 -I- NH 4 SH cloud decks) corresponds to down- 
welling at the deeper, hidden water cloud level. 
Our results suggest that mass flux into the upper 
cloud layers of zones is no t dominated by hori- 
zontal transport, as in the IShowman &: de Pater 
scenario, but is driven by vertical transport from 
below. Jupiter’s circulation in zones therefore 
maintains the same sign of upwelling/downwelling 
across the full 0.5-5 bar weather layer. The same 
sign of upwelling/downwell ing over this large an 
extent was also derived bv Ide Pater et al. ( 2010h 
from 5-micron bright rings around vortices. These 
authors suggested that vortices must extend ver¬ 
tically from at least the 4-7 bar level up to the 
tropopause. 


Large CH 3 D equivalent width requires low wa¬ 
ter cloud opacity. Subsidence of dry upper-level 
air provides the simplest mechanism. Figure 8 
shows that at low latitudes, the clearest deep at¬ 
mosphere (high values of the black curve) occur 
where the 5-/im flux is highest (red curve). The 
very highest flux levels correspond to an atmo¬ 
sphere largely devoid of clouds of all three types, 
H 2 O, NH 3 , and NH 4 SH, again suggesting down- 
welling circulation that spans the full 0.5-5 bar 
range at least. However, the width (in latitude) 
of the 5-^m flux peak is very narrow, compared to 
the width of the NEB/SEB regions with low water 
cloud opacity. Thus, over a broad latitude range 
with little or no water cloud opacity, the atmo¬ 
sphere is characterized by both cloudy and cloud- 
free conditions at the upper levels. Perhaps a two- 
layer circulation model with a vertical flow reversal 
layer between the water and NH 3 / NH 4 SH clouds 
can be relevant on these regional scales, if not for 
the entire planet. 


The rapid “fading” or whitening of the South 
Equatorial B elt between 2008 an d 2010 was in¬ 
terpreted bv lEletcher et al.l ( 201lll as due to en¬ 
hanced upwelling of ammonia-rich air followed by 
condensation. Using CH 3 D line profiles and ab¬ 
sorption lines of NH 3 and H 2 O in the 5-/rm win¬ 
dow, we now have the capability of measuring 
cloud structure and volatile abundances in the 
NEB and SEB at the 4 to 8 bar level. Any fu¬ 


ture changes in the appearance of the belts can 
now be investigated over the full 0.5 to 8 -bar range 
of Jupiter’s troposphere using spectroscopy at 4.6 
and 8.6 /rm. By studying changes in each cloud 
layer separately we will have a much better un¬ 
derstanding of the dynamics below Jupiter’s visi¬ 
ble clouds. 


The high latitudes also show high CH 3 D equiv¬ 
alent widths. Here, more detailed modeling will 
be needed to disentangle geometric effects caused 
by viewing geometry, as well as changes in at¬ 
mospheric scale height in a rapidly-rotating non- 
spherical planet. But the data do suggest a differ¬ 
ent paradigm in deep cloud structure at high lat¬ 
itudes (polewards of 40°). We note, though, that 
5-micron images at high spatial re solution show a 


lot of structure at these latitudes (de Pater et al 


2011 r) . and micro wave images sh ow an overall low 
NH 3 abundance ( de Pater 198^. A combination 
of 5-/rm imaging, microwave imaging, and 5- 
/un spectroscopy of Jupiter’s polar regions would 
be extremely useful to investigate these interesting 
cloud features. 


We now return to the interpretation of 5- 
/un spectra of Jupiter acquired by the KAO, 
Voyager/IRIS, and Galileo/NIMS. It now seems 
clear that the airborne observations of Jupiter 
were flux-weighted by Hot Spots. The abun- 
dances of H 2 O and o ther molecules derived by 
Bioraker et al.l ( 1986al f9l pertain to Hot Spots, 


but not to Jupiter as a whole. The lower bound¬ 
ary for Hot Spots is in fact due to pressure-induced 
H 2 opacity, rat her than an opaque water cloud. 

However, 
does ap¬ 


as proposed by ICarlson et al.l (119921). 


the model proposed by ICarlson et al 


pear to apply to Jupiter’s zones, at least in the 
regions marked in Fig. 8 as candidates for wa¬ 
ter clouds. Similarly, models of Galileo/NIMS 


spectr a of the Equatorial Zone by iDrossart et al 


(|l998r) are probably accurate while m odels of Hot 


Spots that included water clouds (e.g. iNixon et al 
(| 200 lll l will need to be revised. 


This technique provides a constraint on Jupiter’s 
deep water abundance and therefore its 0/H ra¬ 
tio. The base of the water cloud is sensitive to the 
abundance of water because higher abundances 
lead to condensation at deeper levels. The data 
provide the level of the cloud top, not the cloud 
base. Since the top is at higher altitude than the 
base, cloud top constraints provide lower limits to 
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the pressure of the cloud base, or, lower limits to 
the deep abundance of water. A spectrum requir¬ 
ing a water cloud at P > 5 bar would establish a 
supersolar enrichment of water in Jupiter, better 
constraining planetary formation models (Wong 
et al. 2008). 

Figure 4 shows the effect of water cloud pres¬ 
sure level on model fits to a cloudy zone spec¬ 
trum at 32°S. The spectrum is best fit by a water 
cloud with a top between 4-5 bars. The cloud 
base is therefor e found at P > 4 — 5 bar. Fol¬ 
lowing Fig. 1 in Wong etabl ( 2008h . a cloud base 
at 4-5 bar corresponds to 0/H ratios 0.33-1.1 x 
solar (corrected to the new solar 0/H ratio of As- 
plund et al. 2009). Our observations thus pro¬ 
vide a lower limit to Jupiter’s water abundance 
of 0.33-1.1 X solar. This result is consistent with 
the Galileo Probe lower limit of 0.48±0.16x solar 
(Wong et al. 2004), and therefore does not provide 
any new constraint on the water abundance. In fu¬ 
ture work, we will search for spectra that require 
an even deeper water cloud. 


Knowledge of the deep cloud structure permits 
us to retrieve abundances of H 2 O, NH 3 , and other 
5-/rm absorbers more accurately. This will enable 
studies of dynamics below Jupiter’s visible cloud 
layers. Additional observations will permit us to 
observe discrete cloud features such as the Great 
Red Spot and Oval BA. 


The data presented were obtained at the W. 
M. Keck Observatory, which is operated as a sci¬ 
entific partnership among the Galifornia Institute 
of Technology, the University of Galifornia, and 
the National Aeronautics and Space Administra¬ 
tion. The Observatory was made possible by the 
generous financial support of the W. M. Keck 
Foundation. The authors extend special thanks 
to those of Hawaiian ancestry on whose sacred 
mountain we are privileged to be guests. Without 
their generous hospitality, none of the observa¬ 
tions presented would have been possible. We 
also would like to thank Linda Brown for steering 
us to the latest broadening coefficients for GH 3 D. 
This research was supported by the NASA Plan¬ 
etary Astronomy (PAST) Program grant number 
NNX11AJ47G, NNX14AJ43G, NNX15AJ41G, 
and NASA Outer Planets Research Program grant 
number NNX11AM55G. 
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